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Abstract 

We compute the shear viscosity of the unitary Fermi gas above the superfiuid transition 
temperature, using a diagrammatic technique that starts from the exact Kubo formula. 
The formalism obeys a Ward identity associated with scale invariance which guarantees 
that the bulk viscosity vanishes identically. For the shear viscosity, vertex corrections 
and the associated Aslamazov-Larkin contributions are shown to be crucial to reproduce 
the full Boltzmann equation result in the high-temperature, low fugacity limit. The 
frequency dependent shear viscosity r](ui) exhibits a Drude-like transport peak and a 
power-law tail at large frequencies which is proportional to the Tan contact. The weight 
in the transport peak is given by the equilibrium pressure, in agreement with a sum rule 
due to Taylor and Randeria. Near the superfiuid transition the peak width is of the 
order of 0.5 Tp, thus invalidating a quasiparticle description. The ratio 77/s between the 
static shear viscosity and the entropy density exhibits a minimum near the superfiuid 
transition temperature whose value is larger than the string theory bound K/(AirkB) by 
a factor of about seven. 

PACS: 67.10.Jn; 67.85.Lm; ll.30.-j 



1. Introduction 

The remarkable derivation of a simple proportionality between the shear viscosity 77 
and the entropy per volume s in a Af = 4 supersymmetric Yang-Mills theory in the limit 
of infinite 't Hooft coupling A = g 2 N — > 00 by Policastro, Son and Starinets [1] and the 
conjecture by Kovtun, Son and Starinets (KSS) [2] that the ratio 77/s is larger than the 
value h/ (AirkB) found in this limit for all scale invariant, relativistic field theories have 
motivated the search for the 'perfect fluid' which realizes, or at least comes close to, this 
bound [3]. In spite of some theoretical counter-examples [4, 5, 6, 7, 8], the KSS conjecture 
turns out to be valid for all real fluids that are known. In particular, since the velocity of 
light does not appear in the KSS bound, the conjecture is applicable to both relativistic 
and non-relativistic field theories and may even be extended to complex, classical fluids. 
In the case of water, for instance, the minimum value for the ratio 77/s is found close to 
its critical point at 650 K and is only a factor of 25 above the KSS bound [2, 3]. Clearly, H 
is irrelevant for the minimum value of 77/s of water at these temperatures. Yet, as shown 
in Appendix A, there is a simple argument which shows that the viscosity minimum of 
purely classical fluids is not far above that expected from the KSS bound. This begs 
the question what are the necessary conditions for a fluid to be 'perfect' in the sense 
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of a minimum value of rj/s that is limited only by quantum mechanics and, moreover, 
what is the role of scale invariance in this context? As pointed out by Kovtun, Son and 
Starinets [2] and discussed in more detail in a number of recent lecture notes [9, 10, 11] on 
connections between holographic duality and many-body physics, a crucial requirement 
for a fluid to come close to the KSS bound is the fact that it is strongly interacting and 
thus has no well-defined quasiparticles. Indeed, in a situation with proper quasiparticles, 
transport coefficients like the shear viscosity can be computed using kinetic theory. Since 
the lifetime broadening h/r is much less than the average energy fc^T for well-defined 
quasiparticles, the resulting rj/s ^> h/ks is typically far above the KSS bound [3]. By 
contrast, for strongly coupled quantum field theories, the relaxation times are expected 
to be of order h/iksT). Since r\j 1 sT is a characteristic time scale for shear relaxation, this 
immediately implies that rj/s is of order h/ks in the strongly coupled limit. A nontrivial 
example in this context is the standard SU(3) Yang-Mills theory. The associated rj/s ratio 
for the pure gauge theory has been calculated numerically using lattice QCD [12, 13, 14]. 
Its minimum appears near the deconfincment transition temperature and turns out to 
be rather close to the KSS bound. The standard Yang-Mills theory near T c and its 
supersymmctric extension — for which rj/s is independent of temperature and no well- 
defined quasiparticles exist at arbitrary energies — thus have very similar values of rj/s. 
Adding fermions to the pure gauge theory, the ratio rj/s becomes an experimentally 
accessible quantity in high-energy, non-central collisions of heavy nuclei. The observed 
ratio rj/s of the quark-gluon plasma is around 0.4 h/ks, i.e., a factor five above the KSS 
bound [3]. 

In a condensed matter context, generic examples for strongly coupled, finite temper- 
ature field theories for which no quasiparticle description holds are provided by models 
that exhibit a zero temperature critical point [15]. These systems are scale invariant at 
a critical value g c of the coupling g. At finite temperature T, there is a quantum critical 
regime above the critical point, which covers a finite window T > \g— g c \ zu ■ In this regime, 
the thermal energy fc^T is the only energy scale and correlations of the order parame- 
ter exhibit incoherent relaxation with a characteristic time scale = Cft/(fcsT) [11]. 
Here, C is a constant that only depends on the universality class of the quantum phase 
transition. In addition, universal behavior shows up in transport coefficients like the con- 
ductivity and also the shear viscosity in the hydrodynamic regime hkj <C ksT. A concrete 
example is the pseudo-relativistic theory of graphene where the ratio r\j ' s = <& v h/kB has 
recently been calculated within a Boltzmann equation approach. The marginally irrele- 
vant Coulomb interaction in this case gives rise to a logarithmic temperature dependence 
Q n — 0.008 ln 2 (l/T) of the prefactor $^ [16]. This yields a monotonically increasing 
viscosity as the temperature approaches zero within the quantum critical regime. Loga- 
rithmic singularities as T — >• are also present in standard 2d Fermi liquids [17]. 

Here, we consider the shear viscosity for the unitary Fermi gas, a system of attrac- 
tively interacting Fermions at infinite scattering length. The unitary Fermi gas is realized 
experimentally with ultracold atoms in a balanced mixture of, e.g., the two lowest hyper- 
fine levels of 6 Li at a Feshbach resonance and has been studied quite extensively over the 
past few years [18, 19, 20]. It provides an example of a non-relativistic field theory that is 
both scale- and conformally invariant [21]. The underlying quantum critical point in this 
case is the zero density gas at unitarity, as was shown by Nikolic and Sachdev [22]. As 
a consequence of scale invariance, pressure p and energy density e of the gas are related 
by p = 2e/3 [23]. Moreover, the bulk viscosity vanishes at all temperatures [24, 25]. For 
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the shear viscosity, quantitative results so far are only available in the high-temperature 
limit T 3> Tp where a classical description in terms of a Boltzmann equation is possible 
[26], and also deep in the superfluid regime T«T C ~ 0.15 Tp. In the superfluid, a finite 
viscosity arises from phonon-phonon collisions in the normal fluid component, giving rise 
to a rapid increase r)(T) <~ T~ 5 of the viscosity as the temperature approaches zero [27]. 
Since r?(T) ~ T 3 / 2 also increases in the classical limit, both the viscosity and the ratio 
r\l s of the unitary gas necessarily exhibit a minimum, a behavior which is in fact typical 
for any fluid [3]. The major question which will be discussed in this paper, is where 
this minimum appears and what the associated value of the viscosity and entropy is. It 
turns out that the minimum in rj/s for the non-relativistic unitary gas is about a factor 
of seven above the KSS bound, rather close to the value that is found for the relativistic 
quark-gluon plasma. We also show that the minimum in rj implies a lower bound on the 
shear diffusion constant D ni which is about 0.5 h/m. 

In detail, we determine the real part r)(uj) of the frequency-dependent shear viscosity 
from a diagrammatic method that evaluates the stress tensor correlation function in the 
exact Kubo formula. Within a conserving approximation that respects all symmetries and 
the associated conservation laws, we obtain t](lu) in the normal phase and, in particular, 
the static viscosity to entropy density ratio r/(cj = 0)/s. Our results for 77(0;) show a 
Drudc-likc transport peak around lu ~ 0. Its width defines a viscous scattering rate 1/tv, 
which obeys fi/r^ <C ksT at high temperatures, where kinetic theory is applicable. Near 
the transition temperature T c ~ 0.15 Tp to the superfluid phase of the unitary gas, the 
width approaches 0.5 ksTp, thus clearly invalidating a quasiparticle description of viscous 
transport in this regime. The weight W — 2e/3 = p of the Drude peak turns out to be 
equal to the pressure at all temperatures, consistent with a sum rule derived recently by 
Taylor and Randeria [28]. For large frequencies fko > ksT, there is a crossover from the 
Drude peak to an inverse square-root tail ri(uj) <~ 1/y/uJ, whose amplitude is proportional 
to the Tan contact density C [29] . In the high-temperature limit T > Tj? we complement 
our numerical solution of the transport integral equations by an analytical solution to 
leading order in the fugacity z — e 13 ^. We thus confirm that (i) the Kubo formula yields 
exactly the same expression as the Boltzmann equation in this limit, and (ii) the vertex 
corrections, in particular the Aslamazov-Larkin contributions, are crucial even in the 
classical limit and increase the scattering time by a factor of almost three. This resolves 
an inconsistency between the Boltzmann equation result and a previous diagrammatic 
calculation in the high-temperature limit by Bruun and Smith [26]. Finally, we derive 
Ward identities associated with the scale invariance of the unitary gas and show that 
they are obeyed within our approximation for the bulk viscosity vertices at all frequencies 
and temperatures. As a consequence, the bulk viscosity vanishes identically £(cj) = 
[24, 21, 28]. 

The paper is organized as follows: in section 2 we define the microscopic model for the 
unitary gas and discuss its scale invariance. The basic correlation functions that allow 
to calculate the shear viscosity from first principles via a Kubo formula are defined in 
section 3 together with a discussion of hydrodynamic relations and sum rules. In section 
4 we give a qualitative discussion of the temperature dependence of the shear viscosity 
in the classical and in the deep superfluid limit, as well as near the superfluid transition 
temperature. The transport integral equations and their solution within a self-consistent 
Luttinger-Ward formalism are derived in section 5. An exact solution of these equations 
to leading order in the fugacity is presented in section 6, which turns out to reproduce 
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the full solution of the Boltzmann equation in the high temperature limit. In section 7 
we discuss our results for the frequency dependent viscosity in the low temperature — 
but still normal fluid — regime of the unitary gas, where the assumption of independent 
collisions of well defined quasiparticles that underlies the Boltzmann equation breaks 
down. In particular, we calculate the shear viscosity to entropy density ratio rj/s, which 
displays a minimum value slightly above the transition to the superfluid. Conclusions 
and open questions are presented in section 8. There are three appendices: in Appendix 
A, we give an argument based on purely dimensional analysis which explains why the 
inherently quantum mechanical KSS bound applies even to purely classical systems like a 
Lennard- Jones fluid. Appendix B is devoted to the derivation of the Ward identities that 
result from scale invariance and the associated vanishing of the bulk viscosity. Finally, 
in Appendix C we analytically compute the tail coefficient of the shear viscosity at large 
frequencies. 



2. The unitary Fermi gas and scale invariance 



The basic model Hamiltonian that describes a two-component Fermi gas with contact 
interactions of strength g is given by 



H 



C P+9tV- 



qi C P'i C P\ 



(1) 



Here e p = p 2 /2m is the kinetic energy, while <r =t,| denotes the two fermion species 
(hyperfine levels in the context of cold atoms). Since a contact interaction in three 
dimensions needs to be regularized, the coupling strength g(A) is cutoff-dependent. It 
is related to the renormalized physical coupling strength g — 4*7r ft 2 a/m that is fixed by 
the scattering length a via g(A) = g/(l — 2aA/7r). The unitary gas corresponds to the 
special case g = oo of infinite scattering length. In order to see that the system is scale 
invariant at this point, it is convenient to model the contact interaction by an exchange 
of bosons. In a path integral language, the partition function at temperature T = 1//3 
can then be written as a double integral 



= J 



(2) 



over Grassmann fields ?/v, Vv with a =t,| and a complex scalar field -0s, -0s- In the 
case of a Feshbach resonance, this field has a direct physical interpretation as a bound 
fermion pair in the closed channel. The action that describes the interacting Fermi gas 
in this two-channel description is 



S = Jo ^ ' '' '' 



h 2 v 2 



2m 



Vv + i>i 



B 



h 2 v 2 

4m 

+ g($ B ip^i + h 



c. 



(3) 



Here, v is the detuning of the bosonic field, which allows to change the induced 
interaction g ~ —g 2 /v between the fermions. In particular, if g is held at the fixed 
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point corresponding to zero effective range, g = oo is reached at resonance v = 0, where 
the bosonic field is massless. As was shown by Nikolic and Sachdev [22], the action 
S[tp a ,tpB] at this point is invariant under scale transformations x — > Xx and r — > A 2 r. 
The fermionic field ip a (x,T) — > A~ 3 / 2 ?/v(A:c, A 2 t) transforms with the canonical scaling 
dimension, while the bosonic field iPb{x,t) — > A _2 -0s(Ax, A 2 r) acquires an anomalous 
dimension. As a result, the dynamics of the bosons, which arises from the second term 
in equation (3), is irrelevant at unitarity. The description of a finite fermion density 
requires adding a contribution /i ip a ip a to the action. With a fixed value of the chemical 
potential /i, this term explicitly breaks scale invariance. For the unitary gas, however, 
the relation [i ~ £(0)kp between the chemical potential and the Fermi wave vector kp is 
identical to that of an ideal Fermi gas, except for a universal constant £(0) ~ 0.36 — 0.4 
that describes the lowering of the Fermi energy due to the attractive interactions [19, 20]. 
A scale transformation thus leaves the action invariant provided n — > A~ 2 /i is changed 
accordingly to account for the change in density. Alternatively, one may introduce the 
chemical potential through a special choice for the (common) phase of the fields Vv and 
ipB [30]. Scale invariance gives rise to a conserved dilatation current [21] which implies 
the relation e = 3p/2 and the vanishing of the bulk viscosity: during a uniform expansion 
the system remains self-similar and no entropy is produced [25, 24]. In fact, the bulk 
viscosity of the unitary gas vanishes for all frequencies [28] , a result that will be derived 
from an exact Ward identity in Appendix B. A subtle point in this context is the issue 
of whether scale invariance of the action will be present in the full theory, i.e., whether 
the classical symmetry of S survives quantum fluctuations [31, 32]. For the unitary gas 
with equal masses of the two spin components, which is the situation discussed in this 
work, this is believed to be the case. For different masses, however, scale invariance will 
be broken due to the appearance of an Efimov effect, i.e., three- and four-body bound 
states appear even in the absence of a two-body bound state [33, 34]. In particular, 
four-body bound states appear for mass ratios larger than a critical value 13.384 [35]. 
The breaking of scale invariance through the appearance of an Efimov effect has been 
discussed so far mostly for bosons [36] and leads to limit cycle flows in a renormalization 
group treatment [37]. 

3. Kubo formula for the viscosity 

The unitary Fermi gas is in a fluid state at arbitrary temperatures and thus has well 
defined bulk and shear viscosities £ and rj. On a microscopic level, the viscosities are 
the zero frequency limits of linear response functions that may be calculated from first 
principles by the Kubo formula. Specifically, both the shear and bulk viscosities follow 
from the retarded correlation functions of the stress tensor 11^ (in the zero external 
momentum limit), 

x f M (q = Q,uj)=i J dtd 3 a:e <wt fl(t)([II jj -(sB,t),n w (O,0)]) ) (4) 

and its imaginary parts Imx^^u;), which are odd functions of w. In particular, the real 
part of the frequency dependent shear viscosity is determined by the associated positive 
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and even spectral function 



Im y ret (oj) 

Rer](uj) = Xxv ' xyK ' (5) 



and its static limit r/ = lim^^o Rcr;(cj). A completely analogous expression exists for the 
bulk viscosity ((oj) which involves the trace IT, of the stress tensor, 

ReC(w) = ^ . (6) 

The microscopic expression for the stress tensor operator that enters equation (4) for a 
non-relativistic quantum many-body system has been derived by Martin and Schwinger 
[38]. In the case of a scale-invariant short-range potential, where —rd r V(r) = 2V(r), it 
is given by 

j_2 p p 

^i j (q = 0,t) = -J2k i k j c{ (T c ka + I d 3 x I d 3 r r -^2V(r):n t (x)n i (x + r): . (7) 

ka 

The second term explicitly involves the interaction and guarantees that the stress tensor 
satisfies local momentum current conservation as an operator equation (see appendix 
B in Rcf. [21]). In the zero-frequency limit, the kinetic part of the stress tensor (7) 
is the one that appears in the viscous terms of the hydrodynamic equations [3], and 
also as a source term in the Boltzmann equation [39] . In general, the two contributions 
to the stress tensor (7) describe physically different processes for momentum relaxation 
upon insertion into the Kubo formula (4). In particular, the correlation function of two 
interaction contributions describes collisional transport due to interparticle forces, which 
is in fact the dominant contribution in the liquid phase [40]. By contrast, the kinetic 
part is associated with the transfer of transverse momentum due to free particle motion 
and dominates in the gaseous phase. As will be shown in section 5, for the unitary gas 
the interaction part of the stress tensor is important to guarantee that the bulk viscosity 
vanishes, however it gives no contribution to the shear viscosity due to the zero-range 
nature of the interaction. 

Quite generally, momentum current conservation implies that the stress tensor corre- 
lation functions are directly related to the current correlation functions, which are defined 

by 

7${q,w) = i J dtd 3 xe^-^ 6{t) t), j k (0, 0)]) , (8) 

with 

ji(q) = hk * C i- q /2^ C k +q /2,a ( 9 ) 
ka 

the current operator. In the limit q, oj — > 0, these correlation functions acquire a simple 
form, dictated by hydrodynamics. In particular, the existence of a finite shear viscosity 
implies diffusive relaxation of the transverse currents [41, 42]. The transverse part of the 
current correlation spectral function thus has the generic form (for small q, oj) 



2 



OJ oj 2 + (D v q 2 ) 



9 (10) 

2j.m r,2\ 2 



with a diffusion constant D v that is directly proportional to the shear viscosity. Indeed, 
the sum rule 



which quite generally defines the normal fluid density p n [41, 42] immediately implies an 
Einstein relation 



between the shear viscosity and the associated diffusion constant that was first derived by 
Hohenberg and Martin [41]. The shear viscosity can thus be defined both from the stress 
tensor correlation function and, alternatively, from the diffusion constant that appears 
in the transverse current response [42, 28]. For non-superfluid systems, the diffusion 
constant D v is identical to the standard kinematic viscosity v = r]/p. 

For strongly interacting fluids, an exact calculation of correlation functions like the 
ones in equation (4) is hardly possible. It is therefore of considerable interest to find 
constraints that allow to check the validity of approximate calculations. Such constraints 
are derived generically from a short-time expansion of the correlation functions that can 
be expressed in terms of equilibrium expectation values of certain commutators. For 
the viscosities of a Fermi gas with contact interactions, this has been achieved recently 
by Taylor and Randeria [28]. As will be shown in section 5 below, a straightforward 
derivation of the relevant sum rule can be given using a Ward identity due to Polyakov 
[43], who has discussed the behavior of the shear viscosity near the critical point of a 
neutral superfluid within a microscopic approach. This Ward identity, which follows 
from momentum conservation, implies that the frequency-dependent shear viscosity at 
arbitrary values of the scattering length a obeys the sum rule 



Here, in the second form of the equality, we have used the non-trivial relation between 
energy density and pressure away from unitarity, first derived by Tan [44]. The sum 
rule (13) has precisely the form given by Taylor and Randeria [28], however their high- 
frequency tail coefficient C v — C/(10w) is larger by a factor of 3/2 and there is also a 
difference in the prefactor of the contribution proportional to C /a away from unitarity. 
The sum rule implies that r](co) decays like at large frequencies with a prefactor 

that is determined by the Tan contact density C, which is a measure for two fermions 
with opposite spin to be close together [29, 44, 45, 46]. Moreover, the area under the 
frequency-dependent viscosity of the unitary gas is fixed by the equilibrium pressure, 
since the a -1 term is absent at unitarity. Both features are verified with high accuracy 
by our calculations of 7](co) that are detailed in section 7. 

4. Temperature dependence of the shear viscosity 

For the unitary gas, the density n = kp/(3Tr 2 ) fixes both the momentum and energy 
scale, which are the Fermi wavevector Uf and the associated Fermi energy ef = ^bTf = 
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(11) 



V = PnD; 



(12) 




(13) 



h 2 kp I (2m). Since the viscosity has dimensions Hn, purely dimensional arguments require 
the static shear viscosity to be of the form 



rj(T) = hna{6) (14) 

where 6 — T/Tp is the dimensionless temperature scale and a{9) a universal scaling 
function. Remarkably, this function is fixed up to a universal constant 0(1) in the 
high temperature limit. This is based on the counter-intuitive fact that the viscosity 
of a classical gas is independent of its density [47]. Since Tp ~ n 2 / 3 , this requires 
a(6 » 1) ~ 6* 3 / 2 , i.e., a shear viscosity that increases like T 3 / 2 . The same qualitative 
result is obtained from the kinetic theory expression r](T) ~ (p)/a(T) with (p) an average 
momentum and cr(T) a thermally averaged cross section. Since the differential cross 
section for a collision with a given relative momentum p is da/dfl — h 2 /p 2 at unitarity, 
this immediately gives rj ~ p\jT? with px = \f2irmkBT the thermal momentum. Note 
that r\ scales like 1/h 2 even in the classical limit, provided the assumption of zero range 
s-wave scattering remains valid in this regime. 

At low temperatures 8 < 0.15, the unitary Fermi gas is a superfluid. Contrary 
to naive expectations, a superfluid is not a 'perfect fluid' despite the fact that there 
is a vanishing viscosity here. In particular, superfluids do not provide trivial counter- 
examples to the KSS conjecture. Indeed, according to the Landau two-fluid picture, the 
superfluid component has both zero viscosity and zero entropy, so r\j s is undefined at 
T = 0. At any finite temperature, however, a normal component appears, whose entropy 
and viscosity are non-zero. For the bosonic superfluid 4 He this was discussed already by 
Landau and Khalatnikov in 1949 [48]. In the low temperature, phonon-dominated regime, 
they found that the shear viscosity of the normal component grows like T~ 5 because the 
mean free path for phonon-phonon collisions, which are necessary for the relaxation of 
shear, diverges. Specifically, with the assumption that the phonon dispersion has negative 
curvature and thus no Bcliaev decay is possible (this assumption is actually violated in 
superfluid 4 He but is likely to hold for the unitary Fermi gas, see [49]), the calculations 
of Landau and Khalatnikov predict the low temperature shear viscosity to be 

irr, ^ ,^P 2 CI 2 13 (27r) 7 fh Cs \ 9 . . 

Here p n (T) = (2n 2 K/ 45c s )(IcbT / Hc s ) 4 is the normal fluid mass density and u = d\nc s /dlnn 
is the dimensionless strength of the non-linear corrections to the leading-order quantum 
hydrodynamic Hamiltonian which lead to phonon-phonon scattering. The viscosity of 
the normal fluid component thus asymptotically diverges like T~ 5 . As realized by Ru- 
pak and Schafer [27] the same behavior is expected in the unitary Fermi gas. Indeed, 
at temperatures far below T c , the microscopic nature of the superfluid is irrelevant and 
the linearly dispersing Bogoliubov- Anderson phonons are the only excitations that re- 
main. As a result, equation (15) for the viscosity applies also to the unitary Fermi gas, 
provided the exact values of the sound velocity c s and coupling constant u are inserted. 
At unitarity, the sound velocity c s = vp\f(J?> ~ 0.36 vf is directly proportional to the 
Fermi velocity vp with a factor that is determined by the universal Bertsch parameter 
£ ~ 0.36 — 0.4 [19]. Since vf ~ n 1 / 3 as for an ideal Fermi gas, the dimensionless cou- 
pling constant that fixes the strength of the phonon-phonon scattering amplitude has 
the universal value u = 1/3. Together with the standard low-temperature expression 



s = 27r 2 fcs/45 (kBT/hc s ) 3 for the entropy density of a scalar phonon field, the viscosity 
to entropy density ratio 

n£j±°) = A 2 . 15 xio- 5 ^- 8 

s k B 

of the unitary Fermi gas at temperatures far below the supcrfluid transition diverges 
rapidly as (7>/T) 8 . The associated prefactor agrees within 2% with that found from 
a diagrammatic calculation by Rupak and Schafer [27], which is based on the next-to- 
leading terms in the effective field theory of the unitary gas by Son and Wingate [30] (note 
the different prefactor of the entropy in [27]). It is a pure number that only contains the 
Bertsch parameter. Note that due to p n (T) / s(T) — T/c 2 s at very low temperatures, the 
ratio 1)1 ' sT = D n / ' c 2 s is just the characteristic relaxation time t v for shear fluctuations. 
According to equation (15), this relaxation time diverges like t v <~ T~ 9 as T — >■ 0, much 
stronger than the pure thermal energy time scale h/(kBT). In fact, the latter would 
lead to a viscosity to entropy density ratio which approaches a constant at very low 
temperatures, as in the TV = 4 supersymmetric Yang-Mills theory, where no quasiparticles 
exist. 

Near the superfluid transition temperature T c , the universal scaling function a(0) will 
be continuous, yet there will be singularities in higher derivatives. Within the conven- 
tional theory of critical dynamics, the non-analytical behavior of rj = p n D v near T c only 
comes from the thermodynamic singularity in p n (T) [50]. Since p n {T) = p — c' 9(— t)\t\ v 
with d a positive constant and v ss 0.67 the universal critical exponent of the 3D XY- 
model, the conventional theory predicts that a(9) reaches its finite value at T c with an 
infinite slope as the critical temperature is approached from below [51]. A non-analytical 
temperature dependence of the shear viscosity is indeed observed near the A point of su- 
pcrfluid 4 He. In this case, however, a singularity of the form sign(i) \t\" appears on both 
sides of the transition [52] . The associated exponent lo is consistent with the prediction 
co = v of conventional theory below T c , however it has a different value ui w 0.8 above 
the critical temperature. This behavior can be explained within a semi-phenomenological 
approach [53] which is, however, genuine for a liquid state. Thus it cannot be carried 
over to the case of the unitary Fermi gas, for which the precise behavior of a(6) near the 
critical point remains an open problem. 

5. Diagrammatic evaluation of the stress tensor correlation functions 

In order to compute the stress correlation functions (4) we start from the single- 
channel model and derive the corresponding stress tensor in the two-channel model. 
This has the advantage that the potential (fermion interaction) term of the stress tensor 
becomes simply a detuning (mass term) of the bosonic field in the two-channel model. 
We then derive the exact expression (30) for the stress correlation function in the two- 
channel model, including both kinetic and potential terms, in the zero-range limit. This 
expression is evaluated in the self-consistent T- matrix approximation [54, 55]. In the 
two-channel model we compute the fermionic and bosonic self-energies self-consistently 
at the one-loop level but neglect the loop corrections to the Yukawa coupling between 
fcrmions and bosons which appear only at two-loop order (cf. Fig. 1). Variation of these 
coupled single-particle equations with respect to a time-dependent external field results 
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(16) 



in a set of transport equations for the viscosity response functions. A crucial feature 
of this procedure is that it respects the symmetries of the underlying model exactly, 
in particular it obeys scale invariance and thus leads to a vanishing bulk viscosity as 
required. 

Since the continuum model is Galilean invariant, one may expand the stress ten- 
sor operator (7) of the single-channel model (1) into spherical harmonics with angular- 
momentum quantum number I: 

IT = ^2e p y £ (p)4 CT c p(T + 2S m g(A) f d 3 r :n t (r)^(r): (17) 

pa J 

where the spherical harmonics Yi{p) depend on the angle of the vector p = (p, 9, 4>) as 
Ye=o(6, (f>) = 1 and Yi =2 {0, <p) — sin 2 (#) sin(0) cos ((/>). (For convenience we omit the stan- 
dard normalization factor l/\/47r in the spherical harmonics.) For a quadratic dispersion 
the fermion kinetic term admits I = and 1 = 2. Note that for the scale-invariant, 
zero-range model at unitarity the interaction term contributes only for £ = 0! 

In order to compute linear response one has to add a perturbation to the Hamiltonian 
which couples the stress tensor to a time-dependent external field, 



Aff(t)=5^/i/(t)n<(t) = 5^2/i / (t)epy / (p)c+ (T Cp <r + 2ft/ = o(t)ff(A) f d 3 r :n t (r) % (r): . 

J 

(18) 



■pot 



This amounts to the replacement e p M- e p [l + 2he(t)] and g(A) n- <?(A)[1 + 2he{t)Seo] m 
the full Hamiltonian, and likewise in the full action at unitarity 

rP r 

S[h} = / dr [l + 2h t (T)Y e (p)]e p c prT c pa + gip B ipfipi + h.c. 

(19) 



■9 2 7 / 

■tpBWl 



S(A)[1 + 2/i*=o(t) 

The change in the action due to the external field, AS* = S[h] — S, can be parametrized 



as 



AS = 12 \^MX)U„,xx>MX')+i>B{X)U B , xx ,i> B (X')\ (20) 



XX' 

with the coefficient functions 



U a ,xx> = f dTj2 h dr)Ye(p)T^( T XX') , (21) 

For convenience we use a short-hand notation X = (r,r) for the real-space argument 
and K = (k, ie n ) for the Fourier argument, and ^2 X = J d 3 r hr 1 J dr as well as Sxx 1 = 
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S(r — r')M(r — r'). Equations (21), (22) define the bare fermionic and bosonic viscosity 
response vertices at unitarity, 



T%\tX 1 X[) = 



n 2 ViVi 



TO 



5x 1 x' 1 S{t - n) 



si 0) (T*i*D 



(^(AJ-^jq^T-n) for^ = 
|0 for ^ = 2. 



(23) 
(24) 



The response of the grand potential to the external field in terms of the fermionic 
Green's functions G a ,xx> = (Tipa(X)ipl(X')) and bosonic Green's functions Gb,xx' = 
{T^ B {X)^ B {X')) is' 



Sil = tT[(5U a )G a ]+ti[{5U B )G B ] 



(25) 
(26) 



where the trace includes the spin sum and in the second line additionally the i sum. In 
the second line the bosonic Green's function is replaced by the vertex function Txx> = 
—g 2 GB,xx<- Hence, we obtain 



— E, T "'^ xx '> G — + ti + JU)). E s."W> r M . (27, 

In particular, for a static scaling perturbation /i^ = o(t) = h we recover the Tan energy 
formula [29] with the correct UV regularization, 



h=0 



ka 



h 2 k 2 

2m 



«fe CT - = W = 2e - (28) 



where the local limit of the vertex function r 'xx'\x'=x+ nas been expressed in terms of 
the contact density (C4). The Kubo formula (4) in imaginary time can be re-expressed 
in partial-wave components as 



Xi(r) = J d 3 r(Tn,(r,r)n £ (O,0)) 



(29) 



where T denotes time ordering. This stress correlation function can be obtained from 
the second derivative of the grand potential (27), 



s 2 n 



5h e (T)Sh e (0) 



h l= 



tr[T^(0)T^(r)] + tr[s / (0) (0)(&(r) - 4J(r)r(0) 



(30) 



where the last term in the square brackets comes from the explicit hi dependence in 
the second term of equation (27). Note also that this last term is crucial to obtain a 
vanishing bulk viscosity at unitarity a^ 1 = 0, as we will show in Appendix B using 
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the Ward identities that follow from scale invariance. Tg and St are the fermionic and 
bosonic viscosity response functions 



f at (rXX') 
§i(TXX') = 



5G, 



a.XX' 



ST X x 



h t =0 



Sh e (r) 



-g 2 (rn,(r)^ B (X)^(^')>- 



(31) 
(32) 



hi=0 



Note that for the case of the shear viscosity where 5^°' = 0, the viscosity response 

function can be expressed by the L function [54] as fg = LT^°\ and the correlation 

function (30) reduces to the well-known form xt=i(j) = — te[T^ ) =2 (T)LT^ J =2 (0)]. 

The formalism developed so far now allows to derive the sum rule given in equation 
(13), starting from the Ward identity for momentum conservation in the static limit of 
external uj = 0, q — > which reads [43, eq. (5.26)] 



T X y = ~p x dG/dp y . 



(33) 



Here T xy = Ti =2 Yg =2 is the vertex function of the stress tensor (cf. equation (31)). 
Starting from equation (30) for the shear viscosity correlation function, we insert the 
bare shear viscosity vertex T^ 2 (p,ie n ) = h 2 p 2 /m from equation (23) and the full shear 
viscosity vertex T^} 2 {p, ie n ) = —pdG/dp from the Ward identity (33) and obtain at zero 
external Matsubara frequency iuj m = (note that Xxy,xy{iu m ) = X^=2(^ TO )/15 due to 
the angular average of \Yi =2 ] 2 ) 



Xxy,x 



(lLQ r , 



= 0) = 



15^/3 



IK ^ 



pa 

2 x 



5? 



h 2 p 2 dG(p,ie n ) 
P- 



dp 



p 



dn 

pa 

in terms of the fermionic momentum distribution n p = —0 z2 i€ G{p,ic n ). We now 
integrate by parts, with a boundary term at momentum cutoff A, and employ the Tan 
energy formula [29] to express the sum over the momentum distribution by the internal 
energy density e, 



Xxy,xy\i^m 0) — 



5 ^ ' £pl~lp 



h 2 p 5 n p A 
o 



- 1± 

~ ~3~ + 157r 2 m 



2-K 2 m 

pa 

4ft 2 CA h 2 c 



6irma 



(34) 



If the momentum cutoff A is translated into a frequency cutoff f2 = hA 2 /m [28, endnote 
39] we arrive at the sum rule given in equation (13). 

The derivation of the correlation function (30) so far has been completely general and 
exact. The remaining challenge then is to evaluate the viscosity response functions Tg 
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Figure 1: [color online] Diagrammatic contributions to the fermionic and bosonic self-energies. 

and S( within the microscopic model for the unitary Fermi gas. In the following we shall 
do this within the T-matrix approximation. We start from the Dyson equation for the 
fermionic Green's function G^xx' m the- Matsubara formalism, 

G a,xx> = G o^xx> - u «,xx> - ^<x,xx> (35) 

with bare Green's functions Gq{K) = —(ihe n + fi— £*.) and the external field U a .xx' 
from equation (21). In the T-matrix approximation the fermionic self-energy describes 
how fermions scatter off pair fluctuations (cf. Fig. la), 

^a,XX' = G- a .X'X ^XX> ■ (36) 

The Bethe-Salpeter equation for the vertex function which mediates the resonant Fermi- 
Fermi interaction is 

^ xx' — 9{A) 1 — Ub.xx' + G^.xx' G^.xx' ■ (37) 

It contains the inverse bare coupling g(A) _1 = g _1 — mK/ (2ir 2 h 2 ), the external field Ub 
from equation (22), and the bosonic self-energy (cf. Fig. lb). As mentioned above, the 
dynamics of the pair field at unitarity a -1 — > only arises from the excitation of fermion 
pairs while the dynamics of the bosons is irrelevant. Since we consider a balanced gas 
with equal populations of fermion species /j^ = /^j. = fi wc will henceforth drop the spin 
index a . 

The T-matrix approximation for the exact viscosity response vertices is then obtained 
by taking the derivative of the self-consistency equations (35)— (37) with respect to the 
external field hi(r) (cf. Fig. 2); in this way it is guaranteed that the conservation laws 
are satisfied. The amputated viscosity response vertex 

T e (rXX') = J-SffL = £ G -i f e (rYY>) Gy} x , (38) 
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Figure 2: [color online] Diagrammatic contributions to the renormalized viscosity response functions 
Tg and Si. (a) The full fermionic response function Tg has Maki-Thompson and Aslamazov-Larkin 
contributions, in addition to the bare vertex T^ '. (b) The full bosonic response function Sg has two 
contributions which are equivalent for equal fermion masses and populations, in addition to the bare 
vertex . (c) The viscosity correlation function is given by a skeleton diagram with one bare and one 
renormalized viscosity response vertex and full propagators (fermionic and bosonic). In particular, the 
second diagram describes transport via collective modes. 



is given by the variation of the Dyson equation (35) with respect to ht{r), 

T t {rXX') = T^(tXX') + T^ t {tXX') + T e Ah (rXX') (39) 
with the Maki-Thompson and Aslamazov-Larkin vertex corrections 

T™ t (tXX') = T XX , = f e (rX'X) T(XX') , (40) 

T^{tXX') = G x , x = G x , x Sz{rXX') . (41) 

Sht{T) 

The amputated bosonic viscosity response vertex 

S t {rXX') = = J2 T xy Se(rYY>) T Y } X , (42) 

is given by the derivative of the Bethe-Salpeter equation (37), 

S t (rXX') = sf\rXX') -2G XX , T t {jXX') . (43) 

The bosonic viscosity vertex describes the response of the pair field to a scaling or shear 
perturbation. The bare bosonic viscosity vertex serves to remove the UV divergence 
of the particle-particle loop (second term) which is present only for the bulk viscosity 
1 = 0. 
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The self-consistent transport equations for the two-channel model are equivalent to 
those of the standard T-matrix approximation in a purely fermionic description. Specifi- 
cally, the two fermionic vertex corrections T^ MT and T £ AL above are precisely those arising 
in the self-consistent solution of the integral kernel of the L function, as introduced by 
Baym [54, equation (60)]. In physical terms, the first term describes the interaction 
between two fermions by exchange of a single pair while in the second term two pairs are 
exchanged. In the context of calculating the change in conductivity due to superconduct- 
ing fluctuations, these terms are called the Maki-Thompson (MT) and Aslamazov-Larkin 
(AL) contributions respectively [56], a notation that will be used also in the present con- 
text. Note that the second term explicitly includes particle-hole fluctuations, and it is 
also referred to as the "box diagram" in a functional renormalization group approach to 
the thermodynamics of the unitary Fermi gas [57]. These vertex corrections are in fact 
crucial to obtain an approximation which satisfies the conservation laws of the underlying 
model. 

The self-consistent equations (35)-(43) have a structure similar to the equations of 
the Luttinger-Ward approach to the BCS-BEC crossover developed in our previous work 
[58, 59]. In particular, using Fourier transforms and the convolution theorem they be- 
come algebraic equations which afford an efficient numerical solution. The numerical 
calculations are performed in three steps. In the first step the Green's function Gxx' 
and the vertex function V xx' are calculated by solving the self-consistent equations (35)- 

(37) iteratively. Without the external fields U a and Ub the Dyson and Bethe-Salpeter 
equations are diagonal in Fourier space, while the fermionic and bosonic self-energies are 
local in real space. Hence, the coupled equations are solved efficiently by going back and 
forth between real and Fourier space. 

In the second step Gxx 1 and Txx> are used as input for the self-consistent equations 

(38) -(43) to calculate the viscosity response functions TJp, St- Again, the integral equa- 
tions (38) and (42) become algebraic and are solved in Fourier space, while the other 
equations remain local in real space. Note that the spatial Fourier transform between 
radial distances r and radial wavenumber k for the partial-wave component £ is given by 

/>oo 

T t (k) = M-*Y / drr 2 3e (kr)T e (r), (44) 
Jo 

jt poo 

T e (r) = -^j^ dkk 2 ]t (kr)T t (k). (45) 

In the third step the correlation function xti^m) is computed from (30). It is continued 
analytically from the discrete imaginary Matsubara frequencies ioj m to the continuous 
real frequencies u via both the Pade method and a model fit function (cf. section 7). We 
thus obtain the retarded correlation function x^ et (w) = x^(cj) + z x"(cj). Finally, the real 
parts of the viscosities r\ and ( are obtained from the correlation functions for I = 2 and 
1 = according to (cf. equations (5) and (6)) 

Re, H .^M, (46) 

ReCM^^M, (47) 

where the prefactor of 77 comes from the angular integration of the spherical harmonics 
\Xl=2{p)] 2 - Alternatively, one may solve the integral equation directly for real frequencies 
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where the limit ui — > can be taken analytically. In practice, this approach is useful at 
high temperatures, where self-consistency no longer plays a role. 



6. Boltzmann-equation limit 

In the high-temperature limit T ^> Tp the integral equations (38)-(43) can be solved 
by expanding in powers of the fugacity 

z = e^ = t^9- 3/2 + 0(6- 3 ) . (48) 
oV 71 " 

To leading order in z, the pair propagator and on-shell self-energy are given by 

47^3 -3/2 

T ret (k,Q) = -i = + Q(z) (49) 

y/m + 2p- e k /2 

rret, n .8e F erf(y/np/p T ) 

S ( p ,e = e p -») = i — — +0(z). (50) 

In the case of on-shell fcrmions with k = p 1 +p 2 , hVl + 2p = e Pl + £ P2 the pair propagator 
reduces to the well-known scattering amplitude f(q) = i/q at infinite scattering length 
of two particles in vacuum, with relative momentum q. Note that the exact leading- 
order result for the on-shell fermionic self-energy contains a non-trivial error-function 
dependence on the ratio of the momentum p to its thermal value pr that was missing in 
previous studies [60]. It is due to the square-root tail in the pair propagator and gives a 
noticeable correction at thermal momenta p ~ p T . Moreover, this form is indeed crucial 
to fulfill the condition of scale invariance, as will be discussed below. 

The fermionic spectral function in the low fugacity, high temperature limit has most 
of the spectral weight concentrated in the coherent peak at e = e p — p. The peak width 
7 p = ImS ret (p, e) vanishes like efPf/p ~ T -1 / 2 for typical momenta p s=s px, consistent 
with the assumption for the temperature dependence of the relaxation time introduced 
by Bruun and Smith [26]. This implies, in particular, that the fermionic quasiparticles 
become well-defined and thus a Boltzmann equation description is valid in the regime 
<9> 1. 

From a numerical, iterative solution of the integral equations (38)-(43) in the high- 
temperature limit we obtain r\j(hn) = 2.80(1) (T/Tp) 3 / 2 . This fixes the constant in the 
asymptotic behavior a{9) = const 9 3 / 2 at large values of 9 of the universal function intro- 
duced in (14). Within the error bars, the numerical value agrees with that obtained from 
a variational solution of the full Boltzmann equation, using higher Sonine polynomials 
[26, appendix]. The prediction of a simple power-law dependence of the shear viscosity 
r)(T) <~ T 3 / 2 has recently been verified experimentally in a temperature range between 
~ 1.5 and 9 ~ 7 by measuring the expansion dynamics of a unitary gas released from an 
optical trap [61]. Very good agreement has been found also with the expected prefactor, 
thus considerably improving the situation compared to earlier measurements of the shear 
viscosity from the damping of the radial breathing mode [62]. 

Remarkably, the solution of the transport integral equation at high temperatures 
and small frequencies can also be obtained by a completely analytical approach. In 
fact, in the low fugacity limit, one can terminate the iterative procedure after the first 
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Figure 3: [color online] Diagrammatic contributions to the viscosity correlation function Xt(u)) at first 
order in the pair fluctuations: Self-energy (S), Maki-Thompson (MT) and Aslamazov-Larkin (AL) dia- 
grams. 

iteration step (correlation function to first order in the pair propagator) and resum via 
a memory function approach, a method that was developed in the context of electrical 
conductivities by Gotze and Wolfle [63] . The first-order correlation function contains the 
diagrams for self-energy, Maki-Thompson and Aslamazov-Larkin contributions shown in 
Fig. 3. These diagrams are obtained by evaluating the transport equations (38)-(43) 
with the bare viscosity vertices , Sg ^ on the right-hand side, and with bare fermionic 
propagators Go(p, e) with an additional impurity scattering rate 7, which is taken to zero 
after resummation. 

Explicitly, to lowest order both in the fugacity and in the scattering rate, the bosonic 
vertex correction Si in equation (43) is given by 



Here, St in Fourier space (for notation see equation (B.5)) is continued to real frequency 
with one retarded and one advanced boson line, Sg(k,uj + iO, O + u) + iO, O — iO) — > 
Si(k,n + iO,Q — iO) in the limit of vanishing external frequency u) — > 0, to obtain a 
function of wavenumber k and real frequency £1. This bosonic viscosity vertex is used to 
compute the fermionic AL vertex correction in equation (41). 

Next, the fermionic MT and AL viscosity vertices are continued to real frequencies, 
with one retarded and one advanced fermion line, T^(fc,u; + ?0, e + u + iO, e — iO) —> 
Tg(k, e + iO, e — iO) in the limit of vanishing external frequency u — > 0. Finally, since the 
fcrmions are well-defined quasiparticles it is sufficient to evaluate the viscosity vertices 
on-shell (e = e p — /i). Introducing a dimcnsionless momentum variable y — p/y/2rnJtBT — 
\Fk p/pt, the resulting fermionic MT and AL vertex corrections are given by 



From the two fermionic vertex corrections and the self-energy contribution one obtains 
the leading correction to the shear viscosity rj(uj — > 0) according to Fig. 3, 




Se= 2 (k,n + iO,Q-iO) 



32tt7 




(1) (1) (i) _ z 2 {2mf' 2 T 7 / 2 -43 - 3 + 30 _ V2z 2 (2m) 3 / 2 T 7 / 2 

V S +%IT + ^AL- 157r 5/2 n 2 7 2 X lbn^h^ 
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(51) 



which is negative because adding a scattering channel on top of the impurity scattering 
7 lowers the viscosity. Now, the memory function approach effectively assumes that the 
first-order correction r/ 1 ) to the non-interacting viscosity 

(0) _ Z (2 m) ^T^ 
11 ~ 8ttW 7 {b2 > 



is the first term of a geometric series, 



r, = r/°) + V W + ■ ■ ■ « v {0) 



7 ^o (^ 0) ) 2 _ 15gfe| 
rj w ~ 64^/2^ 



e 3 ' 2 . (53) 



Normalized to the particle density, the resulting shear viscosity at high temperatures is 
thus given by 

v = Hn 6 3 ' ' 2 « 2.77 Hn6 3 / 2 . (54) 

64^2 

This analytical result is in perfect agreement with that obtained from the numerical 
solution of the full transport equations above and agrees also with the result obtained 
from the Boltzmann equation, using the standard variation of the distribution function 
Sf oc v x p y [64, 65] . Our evaluation of the Kubo formula within the T- matrix approx- 
imation thus recovers the exact Boltzmann equation result at high temperatures and 
low frequencies. It is important to note that a calculation of the shear viscosity which 
uses a nontrivial Ansatz for the fermionic spectral functions but does not include vertex 
corrections [26], gives the correct T 3 / 2 asymptotic power law at unitarity, however the 
associated prefactor 1.06 is far too small. Indeed, the single-particle scattering time un- 
derestimates the transport scattering time t v by a factor of almost 2.6. The V and AL 
vertex corrections in the Kubo formula are therefore crucial even in the classical limit. 
They are equivalent to solving the full Boltzmann equation, not only its relaxation-time 
approximation. 

Moreover, neglecting vertex corrections violates the Ward identities associated with 
scale invariance and thus leads to a finite bulk viscosity at unitarity. Indeed, it is straight- 
forward to show that only the sum of the three first-order corrections to the bulk viscosity 
(order j~ 2 in the impurity scattering rate) 

C< 1) -<f ) +<a + <ffl-0 (55) 

vanishes identically, as expected for a scale-invariant system. The error-function term in 
the fermionic self-energy (50) is crucial for this cancellation, as well as for satisfying the 
Ward identities derived in Appendix B. 



7. Results for the shear viscosity 



As the temperature is lowered towards the degeneracy temperature Tp, the fugacity 
ceases to be a good expansion parameter, reaching z = 1 at around T rts 0.6 Tp. In a 
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perturbative approach, this regime can be approached from the ideal gas via the virial 
expansion [66, 67, 68]. The known exact values of the second and third virial coefficient 
can then be used for a calibration of measurements of the equation of state of the unitary 
gas [69, 70, 71]. For a calculation of dynamical properties like the viscosity in the regime 
0.15 < 9 < 1, where the gas is degenerate but still in a non-superfluid state, it is neces- 
sary to use a non-perturbative approach that also provides information about dynamical 
correlation functions. As discussed in section 5, this can be achieved by a self-consistent 
Luttinger-Ward formulation of the many-body problem, which allows to determine not 
only thermodynamic properties of the unitary gas at arbitrary temperatures [58] but also 
the fermionic spectral functions [59]. In the relevant regime 0.15 < 6 < 1 just above the 
superfluid transition temperature, they exhibit a substantial broadening, which is of or- 
der of the Fermi energy itself [59]. As a result, the fermionic quasiparticles are no longer 
well-defined. The integral equations (38)— (43) for the viscosity vertices as functions of 
(radial) distance and imaginary time can then only be solved numerically. In practice, 
the self-consistent set of integral equations is iterated until convergence is reached (typ- 
ically after a few steps); then the correlation function Xxy,xy(iu m ) = X£=2(i( j -> m ) / 15 is 
continued analytically to real frequencies. 

As a first step to evaluate the shear viscosity for real frequencies, we perform the 
analytic continuation by the Pade method, using the first 300 Matsubara frequencies 
[59]. The resulting 7](u>) curves are shown in Fig. 4. They exhibit a clear Drude peak, 
which crosses over into a w -1 / 2 tail at large to. Unfortunately, it is difficult to extract a 
precise value for the transport scattering time at large temperatures by this method. To 
obtain reliable results for the frequency-dependent viscosity, we therefore use a simple 
model function of the form 

, W - ^/2h 3 / 2 C r ,m~ 1 / 2 T r] |w m | 3 ^ 2 4H 2 CA , . 

X-(^) = l + TriM + 15^ ' (56) 

It describes the dependence on the Matsubara frequencies found numerically from the 
self-consistent solution of the transport equation very well for arbitrary temperatures 
above T c . The additive contribution in equation (56) which is linear in the momentum 
cutoff ~ CA has been derived in (34). It is, however, irrelevant for the associated viscosity 
spectral function 

LO 1 + (WT,,) 2 -y/mCJ 1 + (cJT r) ) 2 



^ ( Wt v + h 3 / 2 C v T v yJuj/m as lo -> 

1 ft?/ 2 C v / y/muj as lo — > oo 

which exhibits a Drude-like transport peak and a 1/y/uJ tail, as expected on general 
grounds. This viscosity function is in very good agreement with the Pade results shown 
in Fig. 4 and has the advantage that the resulting t v from the fitting function agrees 
perfectly with the exact result in the classical limit. Within this Ansatz, the frequency- 
dependent shear viscosity is characterized by three parameters: the total Drude weight 
W(T), the viscous transport scattering time t v (T) and a tail coefficient C V (T). In prac- 
tice, since both the Drude weight and the tail coefficient are fixed by the equilibrium 
variables p and C via the exact sum rule (13), the transport scattering time remains as 
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Figure 4: [color online] Shear viscosity spectral function r/(uj) for different temperatures (on the left side 
of the plot, temperature decreases from top to bottom; on the right side the order is opposite). The 
analytical continuation is performed by the Pade method. 



the single adjustable parameter. Note that the dc-viscosity in equation (57) has precisely 
the form that is assumed in Maxwell's theory of strongly viscous fluids, with W = Goo 
playing the role of the high-frequency shear modulus [40] . An unexpected feature of our 
model function (57) is the presence of a y/ui singularity at low frequencies, whose weight 
is too small, however, to be seen in Fig. 4. As a result, there is a negative t~ 3 / 2 long-time 
tail in the relaxation of shear stress, similar to that arising in both classical and quantum 
liquids due to mode coupling effects [40, 72]. If we insert the model function into the 
left-hand side of the sum rule (13), 



did 



1 + (WT,) 2 yjvnikl 1 + (WT,) 2 157T-V 



W, 



the second term of the model function and the subtraction of the high-frequency tail 
integrate to zero exactly if C v = C/(15?r), for any value of r^. The frequency integral 
thus yields the Drude weight W, which is equal to the equilibrium pressure p at unitarity 
by the exact sum rule (13). 

As shown in Fig. 5, the total weight of the Drude peak W agrees remarkably well 
with the equilibrium pressure of the unitary gas for all temperatures, from the classical 
limit down to the superfluid transition at T c ~ 0.15 Tp. Here and in the following, 
all thermodynamic properties of the unitary Fermi gas are taken from the equivalent 
Luttinger-Ward calculation of our previous work in Ref. [58]. The viscous transport 
scattering time t v versus temperature is shown in Fig. 6: for large temperatures it 
approaches t v ~ h/sp, consistent with equation (54) and the result obtained from 
kinetic theory [26]. For temperatures near T Cl it is enhanced by a factor of two compared 
with a simple extrapolation of the high temperature result. Qualitatively, there are two 
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Figure 5: [color online] Weight W of the Drude peak in the shear viscosity spectral function ~ 
Wt v /[1 + (ut,) 2 ] as a function of temperature. The pressure p(T) (red line) is from Ref. [58]. 




Figure 6: [color online] Viscous scattering time t v of the Drude peak in the shear viscosity spectral 
function jj(cj) ~ Wt v /[1 + (wt,) 2 ] as a function of temperature. The classical limit (red line) is from 
equation (54). 
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Figure 7: [color online] Static shear viscosity r](w = 0) as a function of temperature. The classical limit 
rj ~ 2.77 (T/T F ) 3 / 2 (red line) is from equation (54). The red diamond indicates T c ~ 0.15 Tp. 



competing effects on the scattering time as the temperature is iowered from the classical 
limit: Pauli blocking reduces the phase space and increases r^, while pairing fluctuations 
enhance scattering and lower t v . In a kinetic theory approximation, both effects nearly 
balance and t v remains essentially at the classical value down to T c [73]. By contrast, 
as is evident from Fig. 6, our evaluation of the Kubo formula within a self-consistent 
T-matrix approximation predicts strong deviations from kinetic theory in the relevant 
regime 6 < 1. Kinetic theory is clearly inapplicable for a degenerate gas and indeed the 
width H/T n w 0.5 ksTp is more than a factor three larger than the thermal energy at the 
lowest temperature T — T c that is studied here. 

The static value r\ = rj{uj = 0) of the shear viscosity is shown in Fig. 7. It exhibits a 
monotonic dependence on temperature, reaching r\ ~ 0.5 tin at the superfluid transition 
temperature, still with a positive slope. As discussed in section 4, one expects a diverging 
positive derivative of r](T) just below T c . As a result, the minimum of the viscosity is 
expected below the superfluid transition, a regime that is not accessible within our present 
approach. The strong rise of the shear viscosity predicted for T <C T c in equation (15), 
however indicates that the minimum value of the shear viscosity for the unitary gas is close 
to the value reached at T c . The dimensionless function a(9) thus has a minimum value of 
order a™ m ~ 0.5, a value that is smaller than any other that is found for non-relativistic 
fluids [3]. Recalling the connection (12) between shear viscosity and the associated 
diffusion constant D v , which is identical with the kinematic shear viscosity v for a normal 
fluid, the minimum value for 77 is equivalent to a minimum value D v ~ 0.5 H/m of the 
shear diffusion constant, which only involves Planck's constant and the particle mass. 
Bounds of a similar form are also expected for other diffusion constants in the unitary 
gas, for instance heat or particle diffusion. In fact, the latter has been measured recently 
by studying the equilibration dynamics after the two spin components are separated in 
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Figure 8: [color online] Tail of the shear viscosity spectral function rj(u) ss h 3 ' 2 C V / \Jmuj as a function 
of temperature. The contact C (red line) is computed numerically according to Ref. [59] and follows the 
high-temperature asymptotics in equation (60). 

a trap [74]. In the degenerate regime, the observed diffusion constant is again of order 
h/m. Concerning heat transport, which so far has neither been calculated nor measured 
for the unitary gas, it is interesting to note that the associated diffusion constant Dt is 
predicted to be identical with the shear diffusion constant D v for non-relativistic fluids 
that exhibit a gravity dual [75]. A comparison of thermal diffusion and shear viscosity 
would thus provide a measure of how close the unitary gas is to such an idealized system, 
whose Prandtl number Pr = D^/Dt is identically equal to one at all temperatures. In 
the non-degenerate regime S > 1 one expects in fact that Pr sw 2/3, a result that is 
essentially universal for dilute classical gases [39] . 

The tail in the viscosity spectral function (57) has the general form 



above a crossover scale u x which we define as the intersection point of the Drude and 
square-root asymptotics. The crossover frequency ui x (T) at large temperatures scales like 
lj x ~ T which agrees with the extension of the Drude peak expected from kinetic theory, 
uj < T [3]. Note that near T c the peak viscosity is still about five times higher than the 
viscosity t](w x ) at the onset of the tail, but the transport peak is nevertheless much less 
pronounced as compared to higher temperatures. In Fig. 8 we plot the tail coefficient 
0^ extracted from our numerical results for the frequency dependent shear viscosity as 
a function of temperature. Apparently, the connection 




H 3 / 2 C, 



+ 0( W " 3 / 2 ) 



(ui > ui x ) 



(58) 



(59) 



between the tail coefficient C V (T) and the Tan contact density C(T) that follows from 
the exact sum rule (13) is fulfilled very accurately at all temperatures (for an alterna- 
tive derivation of the relation (59) see Appendix C). To understand the temperature 
dependence of the tail coefficient, which saturates for low temperatures and decreases as 
<~ 1/T for large temperatures T > Tp, we note that the contact density is defined via 
the asymptotics of the momentum distribution, n a (k) — > C/k 4 for large wavenumbers 
k 3> kp. Alternatively, the contact density can be determined from the vertex function 
(C.4) as h A C = —m 2 r X '=x+ [58]. At high temperatures T ^> Tp this can be evaluated 
analytically and leads to 

_ m _ 4mVr 2 _ 8nW 16k F 

which agrees precisely with the result obtained in Ref. [76] (note the different definition 
of the contact in this work which accounts for an apparent difference by a factor of Air 2 ). 
This asymptotic behavior is in perfect agreement with our numerical results in Fig. 8. 
An alternative way to infer the high frequency behavior of the shear viscosity is based 
on the relation [28] 



3w 3 
q'To Ahq 4 



r){uj) = lim —j lmxpp(q,uj) (61) 



between the frequency-dependent shear viscosity and the mass-density correlation func- 
tion Xpp(?) w )j a relation that is valid at all frequencies. As shown by Son and Thompson 
[77] , the density correlation function at large frequencies 

hfa( ^ 4M)= faW (62) 

is again fully determined by the Tan contact C. The resulting coefficient C v in the high- 
frequency tail of the shear viscosity agrees precisely with that in equation (59) above. 

Our result for the temperature-dependent shear viscosity can now be combined with 
the known value of the entropy density [58] to determine the ratio rj/s in the normal 
fluid regime of the unitary gas. As shown in Fig. 9, this ratio exhibits a very shallow 
minimum around T ss 0.3 — 0.4 Tp, below which rj/s increases very slowly. The precise 
location of the minimum clearly depends sensitively on how accurate the results for 
both the viscosity and entropy are in this regime. On quite general grounds, it is likely 
that the minimum in rj/s is close to the superfluid transition temperature T c ~ 0.15 Tp, 
and that rj/s is monotonically increasing as the temperature is lowered in the superfluid 
regime, eventually crossing over to the steep increase of equation (16) as the temperature 
approaches zero. The fact that our diagrammatic calculation gives results, e.g., for the 
critical temperature and the associated entropy density s ~ 0.7 nks [58] which agree well 
with precise numerical results [78], suggests that our ratio rj/s provides a quantitatively 
reliable estimate, despite the fact that the precise location of the minimum is difficult 
to determine. Granting that our value rj/s ~ 0.6 h/ks for the minimum is close to the 
exact result, we conclude that the ratio rj/s for the unitary Fermi gas remains a factor 
of about seven above the KSS bound, somewhat larger than the experimental estimates 
for the quark-gluon plasma [3] . 
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Figure 9: [color online] Shear viscosity to entropy ratio r)/s (blue circles) in comparison with known 
asymptotes. The dashed red line on the left is the phonon contribution rj/s ~ (T/Tp)~ s in equation 
(16), the solid red line on the right the classical limit (54) divided by the classical entropy from the 
Sackur- Tetrode formula. The red diamond indicates T c ~ 0.15 Tp. 



8. Conclusions 

Using a diagrammatic approach that follows the classic Baym-Kadanoff conserving 
approximation to calculate transport properties of fermionic quantum liquids, we have 
determined the frequency dependent shear viscosity of the unitary Fermi gas in 

its normal phase 9 > 9 C m 0.15. With decreasing temperatures, the Drude peak in the 
viscosity is smeared out and spectral weight is transferred to the tail, precisely 

following the increase in the Tan contact. In the degenerate gas regime 9 < 1, transport 
of transverse momentum cannot be described in terms of a quasiparticle picture because 
the viscosity spectral function is much broader than the thermal energy. Our results for 
r](uj) are in perfect agreement with the exact sum rule (13) which, apart from numerical 
factors, agrees with the one derived by Taylor and Randeria [28]. In particular, the 
weight of the Drude peak, which is effectively the fluid's high-frequency shear modulus, 
is identically equal to the gas pressure at all temperatures above T c , while the coefficient 
C v = C/(157r) of the high-frequency 1/V^ tail agrees with the Tan contact density. 

Our dynamic shear viscosity exhibits a y/oj singularity both at high and low fre- 
quencies. While the high-frequency tail is a consequence of the zero-range interaction, 
providing a direct measure of the contact, the physical interpretation and a detailed the- 
ory of the singularity at low frequencies that is related to the well-known t~ 3 > 2 long-time 
tails due to mode coupling effects remains an open problem. Note that within our for- 
mulation, the sign of the ^fuj singularity is positive, as is found for quantum liquids [72] , 
while the standard long-time tails in classical liquids lead to a negative contribution. 

The proper treatment of symmetries and the inclusion of vertex corrections guarantees 
that our diagrammatic approach leads to a bulk viscosity that vanishes identically. In 
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particular, our approximation exactly satisfies the Tan energy relation (28) and the 
adiabatic sweep theorem [44]. Moreover, it reproduces the high-temperature Boltzmann 
equation result for the shear viscosity by Bruun and Smith [26]. The accuracy of our 
approximation in the most interesting regime just above the supcrfluid transition is 
difficult to estimate a priori. The agreement with the non-perturbative viscosity sum 
rule, however, suggests that our results for the viscosity are quantitatively reliable down 
to the supcrfluid transition. The dc-viscosity is found to decrease monotonically with 
temperature, reaching a value rj ~ 0.5 Hn near T c . Equivalently, the unitary gas near 
T c exhibits a quantum limited shear diffusion constant (or kinematic viscosity) D v ~ 
0.5 h/m. The ratio rj/s between shear viscosity and entropy density exhibits a minimum 
slightly above the superfluid transition temperature. The corresponding value rj/s ~ 
0.6 h/ks is a factor of about seven above the KSS bound, the smallest value found for 
non-relativistic fluids so far. 

There are, of course, a number of problems that remain open. For the unitary gas, 
precise results for the shear viscosity near T c and also immediately below the superfluid 
transition are required to obtain definite predictions for the location and the associated 
value of the minimum in both rj and in rj/s. Similarly, a quantitative understanding 
of current experiments on the inter-diffusion of the two spin components [74] requires 
to study particle transport by a calculation of the conductivity from the longitudinal 
current correlation function. Another transport coefficient that may be determined by 
an extension of our present method is the heat diffusion constant. Its ratio with the 
shear diffusion constant determines the Prandtl number, which is equal to unity in non- 
relativistic, conformal field theories with a holographic dual [75]. 

Finally, it is of interest to study the situation away from unitarity, where scale invari- 
ance is broken and a finite bulk viscosity appears. Since both rj and the ratio rj/s will 
exhibit a minimum as a function of temperature for arbitrary values of l/fc^a along the 
BCS-BEC crossover, this would allow to answer the question whether the minimum found 
at unitarity indeed gives the lowest possible value and thus, whether scale invariance is 
an important ingredient for a minimum in ?j/s in a non-relativistic context. In fact, the 
entropy density at T c has its largest value not at unitarity but near w I [58]. A 

calculation of the viscosity along the complete BCS-BEC crossover has been performed 
in one dimension, where a finite bulk viscosity exists at T = [79] . In this case, the zero 
temperature viscosity is actually smallest in the BEC limit. 

We are grateful for a number of fruitful discussions with Dietrich Belitz, Jens Braun, 
Wolfgang Gotze, John McGreevy, Yusuke Nishida, Mohit Randeria, Achim Rosch, Man- 
fred Salmhofer, Jorg Schmalian, Richard Schmidt, Thomas Schafer and Martin Zwierlein. 
This work has been supported in part by the DFG research unit "Strong Correlations in 
Multiflavor Ultracold Quantum Gases" . 

Appendix A. Viscosity of a Lennard-Jones fluid 

In this Appendix we give an elementary argument, based on dimensional analysis 
only, for the quite surprising observation that the KSS bound which arises from gen- 
uine quantum field theory considerations, even applies to classical fluids, whose viscosity 
minimum apparently cannot depend on K. An example is water, a rather complex fluid, 
that has a minimum viscosity near its critical point at T = 650 K. Since a viscosity has 
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units energy density times seconds, the ratio a n — r)/(hn) is a pure number, which hap- 
pens to be of order one, a v ~ 32 at the viscosity minimum. This appears to indicate 
that quantum mechanics determines the minimum viscosity of water despite the fact 
that the thermal wavelength At near its critical point is much smaller than the average 
interparticle spacing nT 1 / 3 . 

To understand the origin of this surprising fact, consider a simple model for a clas- 
sical fluid, which assumes pairwise central forces that are derived from a Lennard- Jones 
potential 



This potential is, of course, not a realistic model for such complex fluids as water, yet 
it provides a rather good quantitative description of simple fluids like Argon [40]. The 
two parameters a and e of this potential provide the characteristic scales for length and 
energy. As a result, the equation of state p(n,T) = e/a 3 p*(na 3 ,T/e) is completely 
determined by a universal function p*(n*,T*) of the dimensionless density n* = na 3 
and the dimensionless temperature T* = T/e. This provides an explanation for the 
law of corresponding states for different fluids that are reasonably well described by an 
interaction of the Lennard- Jones type. To evaluate dynamic correlation functions like 
the viscosity from a (classical) Kubo formula for this Lennard- Jones fluid requires in 
addition the particle mass m, which fixes a characteristic time scale r = (mo 2 / e) 1 / 2 . By 
purely dimensional arguments the shear viscosity of the Lennard- Jones fluid is then of 
the form 



with a universal function 77* (n*,T*). In fact, equation (A. 2) also applies for the bulk 
viscosity £, which is typically of the same order as 77 in classical fluids [40]. In the 
high-temperature, gaseous limit T* » 1 this function can be calculated from kinetic 
theory, where rj(T) ~ (p) /<x(T) is a ratio between an average momentum and a thermally 
averaged cross section a(T). For large temperatures, the cross section a(T) <~ p T 4 ^™ 
scales with an inverse power of the thermal momentum pt- The exponent is determined 
by the short range repulsive potential $(r — > 0) ~ r~ n . For a Lennard- Jones fluid with 
n = 12, the universal function rj*(n*,T*) <~ (T*) 2 / 3 , therefore, exhibits a power-law 
increase with temperature, independent of density. At low temperatures the viscosity 
will again increase and eventually the response 77 (lo) ~ i/co becomes purely reactive at 
the liquid-to-solid transition, which appears at densities n* > 0.86 [40] (Note that this 
is the density at the freezing point, which is smaller than the melting point density 
due to the first order nature of the transition. Note also that a gas of hard spheres 
with diameter a reaches close packing only at n* = \/2)- Numerical results for the 
dimensionless function rf are available from molecular-dynamics simulations. For the 
Lennard- Jones fluid, they were first performed by Levesque et al. [80]. For example, 
the resulting value is rf (n* = 0.72, T* = 0.84) = 3.29 near the triple point and close 
to the solidification line [81]. The minimum of the viscosity is necessarily reached near 
the gas-liquid transition, which lies between the triple and the critical point. For the 
Lennard-Jones fluid, the dimensionless densities and temperatures at these points are of 
order unity: n* ~ 0.85, T* ~ 0.68 and n* ~ 0.36, T* ~ 1.36, respectively [40]. From 
equation (A. 2), it is thus obvious that the minimum value of the shear viscosity of a 




(Al) 



VL ,( n ,T) = ^- V * („*,T*) 

ST" 



(A.2) 
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Lennard- Jones liquid 

. / 777 f 

Cj n = const ^_ (A.3) 

is fixed by the microscopic ratio y/me/a 2 with a universal prefactor of order unity, that 
unfortunately seems not to have been calculated so far. Clearly, this minimum value does 
not contain H. Indeed, the condition for the applicability of classical physics to describe 
the behavior near the gas-liquid transition of the Lennard- Jones fluid is At <C n -1 / 3 ~ <r, 
which implies that y/me >• hja. The minimum value of the viscosity is therefore much 
larger than h/a 3 ~ tin as long as the standard de Boer parameter 

A* = — %= (A.4) 
cry' me 

is very small. In practice, however, the dc Boer parameter is still appreciable unless one 
considers fluids of very heavy elements like Xc. The associated conversion factor from 
equation (A.3) to a purely quantum limited viscosity n = a v hn is therefore not that 
large. Stated differently, for most fluids the zero-point energy fr? /ma 2 associated with 
the typical length scale of interatomic potentials is not much smaller than the well depth 
e. This is a reflection of the fact that, after all, it is the combination of the quantum 
mechanical van der Waals interaction at large separations and the short range scale a 
that arises from the Pauli principle which determines the ratio s/ms / a . It is, therefore, 
quantum mechanics which eventually fixes the minimum value of the viscosity of classical 
liquids to be of order H times the liquid density. With a corresponding entropy density 
s ~ kgn, this immediately yields a minimum value of rj/s that is not too far above the 
KSS bound, despite the fact that both the viscosity and the entropy may be inferred by 
completely classical considerations. 



Appendix B. Scale invariance and Ward identities 

In this Appendix we show that scale invariance and the associated continuity equation 
for the conserved dilatation current imply Ward identities which relate the fermionic and 
bosonic bulk viscosity vertices to the respective single-particle Green's functions. As 
a consequence the bulk viscosity vanishes identically in our approach as required by 
symmetry. In particular this allows us to check the consistency of the viscosity vertices 
that have been computed explicitly. 

The trace of the stress tensor defined in equations (7) and (17), 

U u {q = Q,t) =U e =o{t) = J2 2s p c U c p°+ I & x I d 3 r2V(r) :n t (x)n^(x + r) : = 2H 

per J J 

(B.l) 

is twice the Hamiltonian in scale- invariant models [82, 21]. The expectation value of this 
operator equation yields immediately 3p = 2e [23], as shown explicitly in equation (28). 
The fermionic bulk viscosity response function Te=o m equation (31) can thus be written 
in imaginary time r as 



T £=0 (rMO = 2(TH{t)^{X 1 )^{X[)) 
28 



(B.2) 



The effects of scale invariance are already incorporated in equation (B.l), and the eval- 
uation of expectation values of the Hamiltonian yields a time derivative [54, equation 
(28)], 

£(H(t)) u = J d\Y\ U ^^G xlx + U x , x ^-G xx ] (B.3) 

where X = (t*,t), X' = (r',r') and the external bilocal held U XX ' couples to fermion 
bilinears via a term ^2 X X' ip'(X) U XX > ip{X') m the action. The coefficient of equation 
(B.3) linear in U X ' lXl is then 

\ ^T^{rX x X' x ) = ± (HirMX^iXD) (B.4) 

d d 
= M(t - t[)—G Xi x[ + fiS{r - T i)^r G Xl x[ ■ 

A Fourier transform of the viscosity response function from imaginary time r to Mat- 
subara frequency iu m and from X, X' to incoming momenta K\ = (k, ie n ) and outgoing 
momenta K[ = (k,ie n + ico m ) yields T^oliUmK^Ki), or 

— — J-e=a(k, tuj m , ie n + iui m , ie n ) 

= (ihe n + ihw m + fi) G(k, ie n + iuj m ) - (ihe n + fi) G(k, ie n ) . (B.5) 

Hence, the bulk viscosity response function is given by a difference of two single-particle 
Green's functions multiplied by their energies. This Ward identity holds exactly and has 
been derived without reference to the T-matrix approximation used above in section 5. 
One can easily check the Ward identity in the non-interacting case where G(k,ie n ) = 
-{ihe n - e k + A*) -1 and f^ (k, iuj m , ie n + iuj m ,ie n ) = G(k,ie n + iuj m )2e k G(k,ie n ). 
As an alternative to starting with equation (B.l), which is already a consequence of 
scale invariance, we have also derived the Ward identity by applying the local conformal 
transformation /3(t) in Rcf. [30] to the generating functional of the Green's functions, with 
the same result. We note that scale plus Galilean invariance imply conformal invariance 
[82] and that both are necessary for a vanishing bulk viscosity [24]. 

A corresponding Ward identity holds for the bosonic viscosity response function (32), 
which is given by the difference of two vertex functions T(K) = T(k, iCl m ), 

= (ihfl m + ihuj m + 2ji) T(fc, i£l m ) — {ifiQ, m + 2/4 ) T(fc, iil m + iui rn ) . (B.6) 

These Ward identities provide a solution to the T-matrix transport equations (38)-(43), 
as one can check by plugging in the Ward identities (B.5) and (B.6) on the right-hand 
side. Thus we have proven that the T-matrix approximation conserves not only charge, 
momentum and energy but also the dilatation current (scale invariance). 

From the Ward identities of scale invariance it follows that the bulk viscosity com- 
puted via equation (30) vanishes for any frequency, temperature or density, 



Re CM = 0. 
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(B.7) 



Indeed, the fermionic Ward identity (B.5) for Ti(K\, K) implies by a shift of the frequency 
that the first term of the stress correlation function (30) vanishes identically. The bosonic 
Ward identity (B.6) implies, again by a shift of frequency, that 

such that also the second term of the correlation function (30) vanishes identically. 

This result has been obtained previously in different ways [25, 24, 21]. For instance, 
it follows in a rather direct manner from a sum rule 

1 I 00 doj C( W ) = — — 2 ( -gZr) (B.9) 

ir J 72irma 2 \dci- 1 J s 

for the bulk viscosity that has been derived by Taylor and Randeria [28]. The fact that 
the derivative dC/d(l/a) is finite at unitarity and that ((co) > is a positive function, 
immediately implies that the bulk viscosity of the unitary gas at a -1 =0 vanishes at 
arbitrary frequencies. 



Appendix C. Relation between contact and tail coefficient 

One can compute the viscosity tail C v defined in equation (58), 



77(^00) = ^=^ (C.l) 



in at least four different ways: from the current correlation function, from the stress 
tensor correlation function, from the relation between the shear viscosity and the density 
correlation function, and from the viscosity sum rule (13). In this Appendix we will 
analytically compute C v at unitarity by the first two methods, thus confirming equation 
(59). The identical result is obtained by considering the high-frequency behavior of the 
density correlation function, as has already been shown in section 7, and from the sum 
rule. 

The trace of the current correlation functions is in general a linear combination of 
shear and bulk viscosities [28] 

xj(q,v) = (\jk(q,u),j k {-q,-u)]) 

= XL{q,u) + 2xt(<?,^) 
'10 



■ 2 

iq 



+ W) • 



At unitarity Q{lo) = 0, and the real shear viscosity is given by 



In the asymptotic high-frequency limit lu — >• 00 only the diagrammatic contributions 
at first order in the pair propagator contribute to the tail. These are the self-energy 
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(S), Maki-Thompson (MT) and Aslamazov-Larkin (AL) diagrams in Fig. 3 (see also [77, 
Fig. 6]), 

X j(q, w) = X S M, w) + xT(Q, w) + X^{q, w) • (C.3) 
In each of these diagrams the vertex function Txx< is replaced by its value at X' = X + , 

-§ 2 (Mb) = Fxx' \x>=x+ = E \ E r ( fc ' in m) = -^r ( C -4) 

which defines the contact density C [29] . The remaining part of the diagram is evaluated 
at zero bosonic wavenumber k = and Matsubara frequency TO = 0, 

r(k,in m )^-^-s(k)S(n m ). (c.5) 

To this order in the pair propagator, the fcrmionic Green's functions are bare Go(p, ie n ) = 
— {ihe n — £ p )~ l and have spectral weight only on-shell. Then to absorb an external 
high-energy perturbation with frequency lo, fermions are excited to very high momenta 
2e p w hui, which are much larger than the Fermi or thermal momentum scales. Hence 
only few-body physics is probed and it is sufficient to evaluate the diagrams in the vacuum 
limit fj, = and set the Fermi functions to zero. 

Specifically, we obtain for the self-energy contribution to the current correlation func- 
tion (C.3) 

a p 

x i E G2 (P' ie n)G{-p, -ie n )G(p + q, ie n + iu m ) + (-q, -iu m ) 

^C.v- P 2 / . \ 

2 V -: , 7TT r + -Q, -iu m ) 

^ Ae p _ (ihw m - e p - s p+ ) 



in 



2 



v 



with the current vertex insertion p squared and shifted momenta p ± = p±q/2. Likewise, 
the Maki-Thompson contribution is 

Xj T (9> lu m) = — 2" EE- 7 ' £ ( P ' :P + ^- : ' fe ( P + 9 ' P ) 
a p 

x i E G (P' ie n)G{-p, -ie n )G(p + q, ie n + iu m )G{-p - q>, -ie n - iw m ) 

= ^ 2E v 2 ( i , i \ 



Expanding to order g 2 and performing the integrals yields the retarded correlation func- 
tion 

rffe, - + ,0) + X r (q, « + i0)= ^^Tjn (C ' 6) 
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and the shear viscosity tail 



h 3 / 2 c 

1s+mt(w -> oo ) = — — == (C.7) 



with tail coefficient C^+mt = C/(Wn). 

In the Aslamazov-Larkin contribution, either one of the two vertex functions is re- 
placed by the contact and the other one is left in place, with the external frequency u 
and momentum q as arguments: 

h 4 C 

Xj L (<?,* w m) = -~^T [Ski=i{q,iu m )] 2 T(q,iuj m ) + (-q,-iu m ) . (C.8) 
The £ = 1 bosonic response vertex with one current vertex insertion pk is given by 
S k ,e=i{q, iu m ) = -2^2jk(p,P + q)^Yl G ( p ' ie n)G(-p, -ie n )G{p + q, ie n + iu m ) 



Pk 



a 2e p _ (ihaj m — e p — e p+ ) 
Integration and analytical continuation to a retarded function yields 



Since [Sk.e=i(q, w)] 2 = —m 3 q 2 /(36TT 2 tvjj) is already of order q 2 it is sufficient to evaluate 
the vertex function (49) at zero momentum, r rot (q = 0,u;) = — 47rift 5 / 2 77i~ 3/ ' 2 cj~ 1 / 2 to 
obtain the AL correlation function 



AL, ^ ■ h3/2C( l 2 



97TTO 1 /2 w 3/2 



and the shear viscosity tail 



The resulting total tail is 



h 3 ' 2 C 

77al(w -> oo) = -— — -=. (C.9) 
30-K\/muj 



n 3 / 2 r 

V (w^oo)= (CIO) 

with tail coefficient C v — C/(15n), in agreement with both the exact sum rule (13) and 
our numerical data. 

An alternative way to compute the tail in the frequency dependent viscosities is to 
consider the stress tensor correlation functions in the limit uj — »■ oo. The calculation is 
analogous to the case of the current correlations, except that there are now stress vertex 
insertions, and one can let the external momentum q — > 0. Note that the interaction 
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contribution to the shear stress operator vanishes for the zero-range model of the unitary 
Fermi gas, cf. section 5. We obtain 



s-iTl,S _ ffll.MT _ C ^n.AL _ n jfi 1 1 ^ 

% -% ^ 30^' % - u y^ AL ) 

in agreement with equation (C.10). For the case of the bulk viscosity, 

c u,s = c umt = C n ,AL = C 

c c 18tt' c 9tt v 7 

such that the total G\ = vanishes in accordance with scale invariance (B.7). 
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